REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1,  REPORT  DATE  (DD-MM-YYYY) 
27-07-2009 


2.  REPORT  TYPE 
Final  Report 


3.  DATES  COVERED  (From  -  To) 

25-Jul-2005  -  31-Dec-2008 


4.  TITLE  AND  SUBTITLE 

Nanoelectronics  of  Carbon  Nanotube  Multi-terminal  Jimctions 

5a.  CONTRACT  NUMBER 

W911NF-05-1-0372 

5b.  GRANT  NUMBER 

5  c.  PROGRAM  ELEMENT  NUMBER 

611102 

6.  AUTHORS 

Madhu  Menon 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f,  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  Kentucky 

University  of  Kentucky  Research  Foundation 

109  KinkeadHall 

Lexington,  KY  40506  -0057 

8.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 

ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

ARO 

U.S.  Army  Research  Office 

P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

1 1.  SPONSOR/MONITOR'S  REPORT 

NUMBER(S) 

48299-EL.2 

12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;DiStributiOn  UnlimitGCt 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

We  have  developed  and  used  a  novel  method  to  investigate  molecular  electronic  device  applications  of  carbon 
nanotube  multi-terminal  junctions  in  the  form  of  T-  and  Y-jimctions.  The  theoretical  methods  included  (i)  large 
scale  simulations  using  accurate  quantum  tight-binding  molecular  dynamics  scheme  and  (ii)  quantum  conductivity 
formalisms.  Both  methods  were  developed  by  the  Pi's  group. 

All  the  structures  investigated  were  fully  relaxed  using  quantum  tight-binding  molecular  dynamics  methods  and 


15.  SUBJECT  TERMS 

Nanoelectronics,  nanotubes,  molecular  dynamics,  quantum  conductivity 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

U 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

U 

17.  LIMITATION  OF 

ABSTRACT 

SAR 


15.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  PAGES  Madhu  Menon _ 

19b.  TELEPHONE  NUMBER 

859-257-8737 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


Nanoelectronics  of  Carbon  Nanotube  Multi- terminal 

Junctions: 

(Final  Report  for  W911NF0510372) 


Madhu  Menon 

Center  for  Computational  Sciences 
University  of  Kentucky 
Lexington,  KY  40506-0045 
Phone:  (859)  257  8737 
Fax:  (859)  323  1029 
E-mail:  madhu@ccs.uky.edu 

July  27,  2009 


Abstract 

Using  the  support  we  have  developed  and  used  a  novel  method  to  investigate  molecular  electronic 
device  applications  of  carbon  nanotube  multi-terminal  junctions  in  the  form  of  T-  and  Y-junctions.  The 
theoretical  methods  included  (i)  large  scale  simulations  using  accurate  quantum  tight-binding  molecular 
dynamics  scheme  and  (ii)  quantum  conductivity  formalisms.  Both  methods  were  developed  by  the  Pi’s 
group.  All  the  structures  investigated  were  fully  relaxed  using  quantum  tight-binding  molecular  dynam¬ 
ics  methods  and  their  current  vs  voltage  (I-V)  characteristics  determined. 


1  Introduction 

Nanoelectronics  of  low  dimensional  materials  provide  an  interesting  study  from  a  fundamental 
physics  of  the  bonding  at  microscopic  level.  These  materials  also  have  useful  potential  in 
technological  applications.  As  semiconductor  device  development  advances  towards  smaller 
devices  and  larger  integration  density,  the  effort  required  for  producing  each  new  generation 
of  memory  and  logic  structures  increases.  Molecular ly  perfect  single- wall  carbon  nanotubes 
(SWCN)  are  the  kind  of  material  possessing  unique  structural  and  electrical  properties  that  are 
expected  to  execute  a  “quantum”  leap  in  the  area  of  nanoscale  electronics. 

The  most  difficult  task  standing  in  the  way  of  the  use  of  SWCN  in  nanoelectric  applications 
is  the  realistic  characterization  of  the  details  of  contact  to  the  outside  terminals  for  the  study 
of  conductance.  The  details  of  metal/nanotube  hetero-j unctions  need  to  be  incorporated  in  the 
overall  circuit  design. 

Theoretical  calculations  can  provide  a  very  valuable  guidance  in  this  area.  With  accu¬ 
rate  determination  of  the  microscopic  structures  of  these  molecular  devices  as  a  first  step, 
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these  calculations  can  subsequently  be  used  to  determine  the  detailed  electronic  properties  of 
these  devices  in  a  manner  useful  to  experimentalists.  Any  realistic  theoretical  determination  of 
structural  properties  of  materials  of  technological  interest  must  involve  the  use  of  molecular  dy¬ 
namics  methods.  Molecular  dynamics  simulations  can  be  used  as  an  “experimental”  technique 
to  accurately  describe  a  specific  material  or  process.  Furthermore,  since  the  goal  is  to  gain  a 
fundamental  understanding  of  structural  and  electronic  properties  at  the  nanoscale  regime,  the 
molecular  dynamics  calculations  have  to  be  based  on  quantum  mechanics. 

Currently,  no  such  simulations  have  been  undertaken.  This  is  due  to  the  difficulty  in  accu¬ 
rately  characterizing  metal-carbon  nanotube  interactions. 

Using  the  ARO  support  we  have  used  quantum  generalized  tight-binding  molecular  dy¬ 
namics  (GTBMD)  method,  developed  by  the  PI,  to  study  nanoelectronic  properties  of  carbon 
nanotubes.  The  theoretical  methods  included  (i)  large  scale  simulations  using  accurate  quantum 
tight-binding  molecular  dynamics  scheme  and  (ii)  quantum  conductivity  formalisms. 

2  Theoretical  Methods 

2.1  Generalized  tight-binding  molecular  dynamics  method 

The  theoretical  simulations  were  performed  using  a  generalized  tight- binding  molecular-dynamics 
(GTBMD)  scheme  developed  by  the  PI  Menon  that  contains  many  state-of-the-art  features[l,  2]. 
Using  this  formalism  we  have  performed  detailed  investigations  of  various  transition  metal  atoms 
interacting  with  single  wall  carbon  nanotubes  (SWGN)  and  find  that  the  curvature  of  the  nan¬ 
otube  and  the  d-band  filling  of  the  d-bands  of  the  transition  metal  atom  has  a  profound  influence 
on  the  bonding  properties  [3].  Using  our  formalism  we  have  shown  the  anomalous  behavior  of 
conductivity  observed  in  experiment  of  SWGN  to  be  due  to  the  effects  of  residual  transition 
metal  catalysts  residing  in  close  contact  with  the  nanotube  walls[4].  Also,  we  have  dynamically 
simulated,  for  the  first  time,  the  action  of  Ni  in  the  carbon  nanotube  growth  process[5]. 

2.2  Quantum  Conductivity 

Although  quantum  conductivity  calculations  of  SWGN  have  been  performed  by  a  number  of 
groups,  use  of  several  simplifying  assumptions  in  the  formalism  may  severely  limit  their  ac¬ 
curacies  while  predicting  useful  device  properties.  A  most  common  practice  is  to  restrict  the 
tight-binding  Hamiltonian  to  only  one  7r-electron  orbital  per  atom.  Furthermore,  most  of  these 
calculations  neglect  the  effects  of  structural  relaxation  altogether.  Even  while  performing  relax¬ 
ation,  most  use  classical  many  body  potentials.  Use  of  tight-binding  Hamiltonian  for  relaxation 
is  very  rare. 

We  have  recently  developed  two  formalisms  for  calculating  the  quantum  conductance  of 
systems  consisting  of  semiconductor  and  transition  metal  systems[6,  7].  They  are  based  on  the 
tight-binding  description  of  the  system.  An  attractive  feature  of  both  formalisms  are  that  the 
tight-binding  Hamiltonian  used  in  the  conductivity  calculations  is  the  same  as  that  used  in 
structural  relaxation.  This  ensures  consistency  in  our  calculations  of  physical  properties. 
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3  Results 

We  have  investigated  the  structural  and  conducting  properties  of  single  wall  carbon  nanotubes 
(SWCNs)  in  embedded  as  well  as  in  side  contact  with  metal  leads  is  obtained  using  efficient 
formalisms  based  on  spin  polarized  tight-binding  formulation  incorporating  full  consideration  of 
s,  p  and  d  basis  sets  for  carbon  and  metal  atoms.  The  full  structural  relaxation  of  the  combined 
SWCN  and  metal  system  is  found  to  be  essential  for  realistic  characterization  of  conductivity. 
Metal  contacting  the  nanotube  and  forming  metal  layer-SWCN-metal  layer  (ML-SWCN-ML) 
system  consists  of  several  layers,  ideally  infinite.  A  few  of  the  metal  layers  are  included  as  part  of 
the  “extended  molecule” .  These  consist  of:  (i)  Metal  layers  in  direct  contact  with  the  nanotube 
and,  (ii)  a  few  layers  that  are  neighboring  to  the  SWCN.  The  “extended  molecule”  consists  of 
nanotube  and  these  two  sets  of  metal  layers.  The  remaining  metal  layers  form  the  two  metal 
leads  (MLs)  each  of  which  is  semi-infinite.  In  the  present  applications,  the  extended  molecule 
consist  of  800-1100  atoms.  The  entire  extended  molecule  is  allowed  to  fully  relax  through  MD 
simulations  with  no  symmetry  constraints. 

The  MLs  are  assumed  to  be  of  bulk  Ni  metal  in  the  <  001  >  orientation  relative  to  the 
nanotube  axis  (taken  to  be  the  z-axis)  in  the  embedded-end  contact  configuration,  while  in  the 
side-contact  configuration  the  Ni<  001  >  direction  is  perpendicular  to  tube  axis. 

In  the  embedded  end-contact  configuration,  the  open  ends  of  the  SWCN  are  buried  into  the 
ML-surface  in  a  way  that  ML-atoms  may  bond  from  both  the  exterior  as  well  as  the  interior 
part  of  the  tube  ends.  The  relaxation  results  in  the  distortions  of  atomic  positions  for  both  Ni 
and  C  atoms  in  the  contact  region.  To  better  illustrate  this,  Ni  atoms  from  the  top  panel  in  Fig. 
1  have  been  removed  to  obtain  the  nanotube  in  the  middle  panel  that  shows  the  distortions  of 
C  atoms  at  the  ends.  The  bottom  panel  shows  distortion  free  nanotube  for  comparison. 

The  calculated  I-V  curves  for  three  embedded-contact  configurations  along  with  the  Ni 
removed  SWCN  system  (Fig.  1,  middle  panel)  are  shown  in  Fig.  2.  These  were  obtained 
by  assuming  symmetric  biasing  i.e.,  VL=-Vi{=V.  The  corresponding  transmission  functions, 
Tcr=Mp(E),  are  plotted  in  the  inset.  From  this  it  can  be  seen  that,  for  (10,0)  SWCN  embedded 
in  MLs,  noticeable  changes  are  induced  in  To-=up(E)  in  going  from  3  to  5  metal  layers.  The 
changes  refer  to  the  location,  the  width  and  the  strength  of  the  resonance  peaks  of  To-=up(E); 
they  are  especially  pronounced  for  energies  below  the  Eermi  energy  E^t’  (set  to  zero),  where  the 
To-=Mp(E)  for  the  system  with  the  3  metal  layers  appears  larger  than  that  of  the  5  metal  layers 
and  exhibits  a  transmission  gap  at  higher  energies  when  compared  to  those  found  in  the  case 
of  the  5  metal  layer  contact.  By  further  increasing  the  thickness  of  the  metal  layers  which  is 
incorporated  in  the  extended  molecule,  i.e.,  in  going  from  5  metal  layers  to  7  metal  layers,  it  can 
be  seen  from  Eig.  2  that  the  differences  in  the  transmission  features  of  the  5  and  the  7  metal 
layer  systems  are  not  significant,  indicating  that  convergence  has  been  attained  in  To-=np(E) 
with  respect  to  the  contact  width. 

We  find  that  convergence  with  respect  to  the  number  of  the  metal  lead  (ML)  atoms  in  contact 
with  the  SWCN  is  found  to  be  very  critical.  Our  results  indicate  that,  in  order  to  maximize 
device  efficiency,  one  needs  to  use  ML-SWCN  systems  with  a  minimal  ML-SWCN  contact-width 
to  SWCN-length  ratio.  If  this  ratio  is  large  enough,  the  SWCN  cannot  be  seen  independently 
of  its  contacts  and  the  ML-SWCN-ML  system  behaves  as  a  molecular  wire.  Additionally,  the 
ML-SWCN-ML  complex  is  found  to  have  spin-selective  transport  properties  for  certain  bias 
range,  making  it  a  promising  candidate  for  use  as  a  spin  valve  in  nanoelectronic  devices. 

We  have  performed  detailed  quantum  conductivity  calculations  for  complex  multi-terminal 
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Figure  1:  (Color  online)  Fully  relaxed  geometry  of  a  (10,0)  SWCN  in  embedded  end-contact 
configuration  in  Ni  (top  panel)  used  in  I-V  calculations.  The  middle  panel  is  obtained  from 
removing  all  Ni  atoms  from  the  top  panel.  The  bottom  panel  shows  a  fully  relaxed  isolated 
(10,0)  SWCN. 
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Figure  2:  (Color  online)  Transmission  functions  for  spin-up  electrons  (upper  inset)  in  units  of 
{e^jh)  and  I-V  characteristics  (current  in  fiA)  for  symmetric  bias  configurations  for  a  (10,0) 
SWCN  in  three  embedded  end-contact  configurations:  (1)  In  Ni-leads  consisting  of  7  (red), 
5  (blue)  and  3  (green)  Ni  layers.  For  comparison,  the  corresponding  results  for  the  same 
nanotube  after  removing  all  Ni  atoms  (Fig.  1,  middle  panel)  are  also  shown  (black  curve).  The 
corresponding  transmission  coefficients  are  shown  by  the  solid  curves  in  the  inset. 
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Figure  3:  Fully  relaxed  geometries  of  two  soldered  multi-branched  junctions  made  of  three 
SWCN-pieces  soldered  in  the  middle  by  a  Ni  cluster.  The  soldering  in  the  right  leads  to  a 
(10,10)-(18,0)-(10,10)  SWCN  T-junction  while  the  soldering  in  the  left  leads  to  an  (8,0)-(8,8)- 
(8,0)  SWCN  Y-junction. 


SWCN  junctions  created  by  nano-” soldering”  with  transition  metal  atoms.  In  particular,  the 
three-terminal  junctions  are  useful  in  that  they  can  be  used  in  device  applications  in  which 
the  3rd  terminal  can  be  used  for  controlling  the  switching  mechanism,  power  gain,  or  other 
transisting  applications  that  are  needed  in  any  extended  molecular  electronic  circuit. 

In  Fig.  3  we  show  fully  relaxed  geometries  of  a  (8,0)-(8,8)-(8,0)  SWCN  Y-junction  and 
a  (10,10)-(18,0)-(10,10)  SWCN  T-junction  created  by  soldering  three  SWCN  branches  by  Ni 
clusters.  All  outer  ends  of  each  of  the  three  arms  of  the  T-  and  Y-junctions  are  taken  to  be  in 
contact  with  the  transition  metal  leads  consisting  of  Ni  atoms  (not  shown). 

The  I-V  curves  obtained  for  the  Y-junction  in  Fig.  3  is  shown  in  Fig.  4.  The  current 
directions  and  the  voltages  on  the  three  arms  are  indicated  in  in  the  upper  inset.  The  main 
graph  shows  the  calculated  currents  in  the  three  arms  of  the  soldered  Y-junction  as  a  function  of 
the  bias  voltage  V3.  The  voltage  configuration  for  this  plot  has  been  set  to  Vi=V2=0.0  V.  This 
setup  makes  the  Y-junction  a  two  terminal  device  for  the  investigation  of  rectifying  behavior. 
As  seen  in  the  figure,  there  is  a  noticible  increase  in  the  current  for  positive  values  of  the  bias 
voltage  V3,  while  for  negative  values  of  V3  the  current  is  negligible.  The  I-V  characteristics, 
thus,  display  a  distinct  asymmetry  and  rectifying  behavior.  While  this  behavior  is  similar  to  the 
non- soldered  case,  although  with  significantly  lower  values  of  the  current  (lower  inset),  a  striking 
difference  can  be  observed.  The  rectification  direction  for  the  soldered  Y-junction  appears  to 
undergo  a  reversal  about  V3=0. 

While  our  results  indicate  that  soldering  of  multi-branched  SWCNs  in  the  form  of  T-  or 
Y-SWCNs  fails  to  reproduce  the  transport  properties  of  the  corresponding  unsoldered  junctions 
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Figure  4:  The  I-V  curves  obtained  for  the  soldered  Y-junction  in  Fig.  3  (left).  The  current 
directions  and  the  voltages  on  the  three  arms  are  indicated  in  in  the  upper  inset.  The  main 
graph  shows  the  calculated  currents  in  the  three  arms  of  the  soldered  Y-junction  as  a  function 
of  the  bias  voltage  V3.  The  voltage  configuration  for  this  plot  has  been  set  to  Vi=V2=0.0  V. 
The  I-V  characteristics  show  asymmetry  and  rectification.  The  corresponding  I-V  curves  for 
the  non-soldered  Y-junction  is  shown  in  the  lower  inset. 
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in  general,  they,  however,  confirm  another  important  feature:  Namely,  they  demonstrate  that 
the  rectifying  and  switching  properties  of  the  Y-  and  T-junctions  are  neither  the  result  only 
of  the  bias  configuration  nor  the  characteristics  of  the  contacts  between  the  tubes  and  the 
metallic  leads.  Because  if  they  were  so,  soldering  could  not  have  such  a  dramatic  effect  on  the 
rectification  properties.  The  present  results  demonstrate  that  the  transport  properties  of  the  Y- 
and  T-junctions  are  strongly  depended  on  their  inherent  properties  as  previously  demonstrated 
in  our  earlier  works  for  the  unsoldered  junctions. 

This  inherent  characteristic  depends  on  the  symmetry  of  the  “spacer”  region  (i.e.,  the  region 
that  joins  the  three  branches)  and  the  way  this  is  bonded  to  the  branches.  In  the  absense  of 
this  inherent  characteristic,  the  spacer  of  the  soldered  junctions  reduces  the  transport  pathways 
which  deviate  from  the  straight  trajectories. 

Our  results  have  demonstrated  that  soldering  reduces  the  current  as  well  as  the  Ohmic 
character.  More  importantly,  for  multi-terminal  SWCN  junctions,  soldering  can  significantly 
alter  the  rectification  properties,  even  reversing  its  direction  due  to  the  destruction  of  the 
symmetry  of  the  junction’s  spacer  region. 
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